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Abstract 

Hamiltonian light-front quantum field theory constitutes a framework for the non-perturbative 
solution of invariant masses and correlated parton amplitudes of self-bound systems. By choos- 
ing the light-front gauge and adopting a basis function representation, we obtain a large, sparse, 
Hamiltonian matrix for mass eigenstates of gauge theories that is solvable by adapting the ab initio 
no-core methods of nuclear many-body theory. Full covariance is recovered in the continuum limit, 
the infinite matrix limit. There is considerable freedom in the choice of the orthonormal and com- 
plete set of basis functions with convenience and convergence rates providing key considerations. 
Here, we use a two-dimensional harmonic oscillator basis for transverse modes that corresponds 
with eigensolutions of the soft- wall AdS/QCD model obtained from light-front holography. We 
outline our approach and present illustrative features of some non-interacting systems in a cavity. 
We illustrate the first steps towards solving QED by obtaining the mass eigenstates of an electron 
in a cavity in small basis spaces and discuss the computational challenges. 

PACS numbers: 11.10.Ef,11.15.Tk 



I. INTRODUCTION 



Non-perturbative Hamiltonian light-front quantum field theory presents opportunities 
and challenges that bridge particle physics and nuclear physics. Major goals include pre- 
dicting both the masses and transitions rates of the hadrons and their structures as seen in 
high-momentum transfer experiments. Current focii of intense experimental and theoretical 
research that could benefit from insights derived within this Hamiltonian approach include 
the spin structure of the proton, the neutron electromagnetic form factor, the generalized 
parton distributions of the baryons, etc. 

Hamiltonian light-front field theory in a discretized momentum basis hj and in transverse 
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lattice approaches [2, y, |J] have shown significant promise. We present here a basis-function 
approach that exploits recent advances in solving the non-relativistic strongly interacting 
nuclear many-body problem . We note that both light-front field theory and nuclear 
many-body theory face common issues within the Hamiltonian approach - i.e. how to (1) 
define the Hamiltonian; (2) renormalize to a finite space; (3) solve for non-perturbative 
observables while preserving as many symmetries as possible; and, (4) take the continuum 
limit. In spite of the technical hurdles. Ken Wilson has assessed the advantages of adopting 
advances in quantum many-body theory and has long advocated adoption of basis function 
methods as an alternative to the lattice gauge approach 

There are three main advantages of Hamiltonian light-front quantum field theory moti- 
vating our efforts to overcome the technical hurdles. First, one evaluates experimental ob- 
servables that are non-perturbative and relativistically invariant quantities such as masses, 
form factors, structure functions, etc. Second, one evaluates these quantities in Minkowski 
space and, third, there is no fermion doubling problem. 

We begin with a brief overview of recent advances in solving light nuclei with realistic 
nucleon-nucleon (NN) and three-nucleon (NNN) interactions using ab initio no-core methods 
in a basis function representation. Then, we introduce our basis function approach to light- 
front QCD within the light-front gauge. Renormalization/regularization issues are also 
addressed. We present illustrative features of our approach with the example of cavity- 
mode QED and sketch its extension to cavity-mode QCD. For a specific QED example, we 
work in small basis spaces and solve for the mass eigenstates of an electron coupled to a 
single photon in a transverse harmonic oscillator cavity. 



The present work is an expanded version of a recent paper where we provided an initial 
introduction to our approach j^. 



II. NO CORE SHELL MODEL (NCSM) AND NO CORE FULL CONFIGURA- 
TION (NCFC) METHODS 

To solve for the properties of nuclei, self-bound strongly interacting systems, with realistic 
Hamiltonians, one faces immense analytical and computational challenges. Recently, ab 
initio methods have been developed that preserve all the underlying symmetries and converge 
to the exact result. The basis function approach that we adopt here [s], [g] is one of several 
methods shown to be successful. The primary advantages are its flexibility for choosing 
the Hamiltonian, the method of renormalization/regularization and the basis space. These 
advantages impel us to adopt the basis function approach in light-front quantum field theory. 
While non-relativistic applications in finite nuclei restrict the basis to a fixed number of 
fermions, we introduce here the extension to a flexible number of fermions, antifermions and 
bosons. 

Refs. jsj and provide examples of the recent advances in the ab initio NCSM and 
NCFC, respectively. The former adopts a finite basis-space renormalization method and 
applies it to realistic nucleon-nucleon (NN) and three-nucleon (NNN) interactions (derived 
from chiral effective field theory) to solve nuclei with Atomic Numbers A = 10 — 13 j^. 
Experimental binding energies, spectra, electromagnetic moments and transition rates are 
well-reproduced. The latter adopts a realistic NN interaction that is sufficiently soft that 
renormalization is not necessary and binding energies obtained from a sequence of finite 
matrix solutions may be extrapolated to the infinite matrix limit. Owing to uniform con- 
vergence and the variational principle, one is also able to assess the theoretical uncertainties 
in the extrapolated result. One again obtains good agreement with experiment. 

It is important to note the analytical and technical advances made to solve these problems. 
First, non-perturbative renormalization has been developed to accompany these basis-space 
methods that preserve all the symmetries of the underlying Hamiltonian including highly pre- 
cise treatments of the center-of-mass motion. Several schemes have emerged with impressive 
successes and current research focuses on detailed understanding of the scheme-dependence 



of convergence rates (different observables converge at different rates) 10|. Second, large 



scale calculations are performed on leadership-class parallel computers, at Argonne National 
Laboratory and at Oak Ridge National Laboratory, to solve for the low-lying eigenstates 
and eigenvectors as well as to carry out evaluation of a suite of experimental observables. 
For example, one can now obtain the low-lying solutions for A = 14 systems with matrices 
of dimension one to three billion on 8000 to 50000 processors within a few hours of wall- 



Ill and the computers are growing 



clock time. Since the techniques are evolving rapidly 
dramatically, much larger matrices are within reach. 

In a NCSM or NCFC application, one adopts a 3-D harmonic oscillator for all the par- 
ticles in the nucleus (with harmonic oscillator energy Q), treats the neutrons and protons 
independently, and generates a many-fermion basis space that includes the lowest oscillator 
configurations as well as all those generated by allowing up to N^ax oscillator quanta of ex- 
citations. The single-particle states are formed by coupling the orbital angular momentum 
to the spin forming the total angular momentum j and total angular momentum projection 
rrij. The many-fermion basis consists of states where particles occupy the allowed orbits 
subject to the additional constraint that the total angular momentum projection Mj is a 
pre-selected value. This is referred to as the m-scheme basis and, in a single run, one ob- 
tains eigenstates with total angular momentum J > Mj. For the NCSM one also selects a 
renormalization scheme linked to the many-body basis space truncation while in the NCFC 
the renormalization is either absent or of a type that retains the infinite matrix problem. In 
the NCFC case jsl, one extrapolates to the continuum limit as we now illustrate. 

We show in Fig. [1] results for the grou nd state of as a function of Nmax obtained 



with a realistic NN interaction, JISP16 



12|. The smooth curves portray fits that achieve the 



desired independence of N^ax and Vt so as to yield the extrapolated ground state energy. Our 
assessed uncertainty in the extrapolant is about 2 MeV and there is rather good agreement 
with experiment within that uncertainty. The largest cases presented in Fig. [T] correspond 
to Nmax = 8, where the matrix reaches a basis dimension near 600 million. Nmax = 10 
produces a matrix near 8 billion and its lowest eigenvalues have now been solved at two 
values of Q. These Nmax = 10 results follow closely the curves shown in Fig. [1] and will be 
presented elsewhere. 




FIG. 1: (color online) Calculated ground state energy of ^^C for N^ax = 2—8 (discrete points) at 
selected values of the oscillator energy, Vt. For each Vt, the results are fit to an exponential plus 
a constant, the asymptote, which is constrained to be the same for each curve[6|]. We display the 
experimental ground state energy and the common asymptote. 



III. CHOICE OF REPRESENTATION FOR LIGHT FRONT HAMILTONIANS 



It has long been known that light-front Hamiltonian quantum field theory has similarities 
with non-relativistic quantum many-body theory. We further exploit this connection, in 
what we will term a "Basis Light Front Quantized (BLFQ)" approach, by adopting a light- 
front single-particle basis space consisting of the 2-D harmonic oscillator for the transverse 
modes (radial coordinate p and polar angle 0) and a discretized momentum space basis for 
the longitudinal modes. Adoption of this basis is also consistent with recent developments 
in AdS/CFT correspondence with QCD Is], 14, 15, 16, 17]. In the present application 



to the non-interacting problem, we will adopt periodic boundary conditions (PBC) for the 
longitudinal modes and we omit the zero mode. For our light-front coordinates, we define 

^) and coordinate pair (p, 0) are the usual cylindrical coordinates 
in [x^^x"^). The variable x'^ is light-front time and x~ is the longitudinal coordinate. We 
adopt x"*" = 0, the "null plane" , for our quantization surface. 



The 2-D oscillator states are characterized by their principal quantum number n, orbital 
quantum number m and harmonic oscillator energy Q. It is also convenient to interpret the 
2-D oscillator as a function of the dimensionless radial variable a/Mq^p where Mq has units 
of mass and p is the conventional radial variable in units of length. Thus, the length scale 
for transverse modes is set by the chosen value of \/MqQ. 

The properly orthonormalized wavefunctions, ^n,m{p,(p) = {p<t>\nm) = fn,m{p)Xm{4>)y 
are given in terms of the Generalized Laguerre Polynomials, L^^iMonp"^), by 

= ^e^"^^ (2) 
with eigenvalues En,m = (2ri + |m| + 1)^2. The orthonormalization is fixed by 

/•oo /*27r 

{nm\n'm') = / pdp dcp ^n,m{p,4>y ^n',m'{P:4>) = ^n,n' ^m,m' (3) 

which allows for an arbitrary phase factor e*" that we have taken to be unity. One of 
the significant advantages of the 2-D oscillator basis is the relative ease for transforming 
results between coordinate space and momentum space. That is, the Fourier transformed 
wavefunctions have the same analytic structure in both coordinate and momentum space, a 
feature reminiscent of a plane-wave basis. 

In order to gain an impression of the transverse modes in our light-front basis, we present 
in Figs. [2] and [3] snapshots of selected low- lying modes. As one increases the orbital quantum 
number m, pairs of maxima and minima populate the angular dependence of the basis 
function. Also, as one increases the principal quantum number n, additional radial nodes 
appear as evident in the progression from Fig. [2] to Fig. [31 

To provide a perspective on the full 3-D basis, we introduce longitudinal modes ipj defined 
on —L < x~ < L with both periodic boundary conditions (PBC) and antiperiodic boundary 
conditions (APBC). We also introduce purely real form to be used in a figure below: 
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FIG. 2: (color online) Modes for n = of the 2-D harmonic oscillator selected for the transverse 
basis functions. The orbital quantum number m progresses across the rows by integer steps from 
in the upper left to 4 in the lower right and counts the pairs of angular lobes. Amplitudes as 
well as X-axis and y-axis coordinates are in dimensionless units. 

where k = 1,2,3,... for PBC (neglecting the zero mode) and k = |,|,|,... in Eq. |4] for 
APBC. Similarly, k = 1,2,3 in Eq. [5] for reflection antisymmetric amplitudes with box 
boundary conditions (amplitude vanishes at x~ = ±L). The full 3-D single particle basis 
state is defined by the product form 

For a first illustration, we select a transverse mode with n = l,m = joined together 
with the k = ^ longitudinal APBC mode of Eq. H] and display slices of the real part of this 
3-D basis function at selected longitudinal coordinates, x~ in Fig. HI For comparison, we 
present a second example with Eq. [5] for the longitudinal mode in Fig. [51 Our purpose in 
presenting both Figs. H] and E] is to suggest the richness, flexibility and economy of texture 
available for solutions in a basis function approach. Note that the choice of basis functions is 



FIG. 3: (color online) Modes for n = 1 of the 2-D harmonic oscillator selected for the transverse 
basis functions. The orbital quantum number m progresses across the rows by integer steps from 
in the upper left to 4 in the lower right and counts the pairs of angular lobes. Amplitudes as 
well as X-axis and y-axis coordinates are in dimensionless units. 

rather arbitrary, including which boundary conditions are imposed, except for the standard 
conditions of orthonormality and completeness within the selected symmetries. 

Although our choice of basis functions is not dictated by theory, it is buttressed by 
the phenomenological success of the "soft- wall" AdS/QCD model jisl . [l^ which uses a 
harmonic oscillator potential in the fifth dimension of Anti-de Sitter space to simulate color 



confinement. As shown in ref. 
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ography" [l^ to transform 



17| to a corresponding bound- 



15| one can use "light-front 
the bound state equations for the wavefunction in AdS space 
state equation in physical space at fixed light-front time r. The resulting light-front equation 
is similar in form to the Schrodinger radial wave equation at fixed t which describes the 
quantum-mechanical structure of atomic systems. However, the formalism at fixed light- 
front time is relativistic and frame independent. Thus, for the specific example of a qq pair, 
one obtains a relativistic wave equation applicable to hadron physics, where the light-front 
coordinate ( = b±y^x{l — x) plays the role of the radial variable r of the nonrelativistic 
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FIG. 4: (color online) Transverse sections of the real part of a 3-D basis function involving a 
2-D harmonic oscillator and a longitudinal mode of Eq. [5] with antiperiodic boundary conditions 
(APBC). The quantum numbers for this basis function are given in the caption. The basis function 
is shown for the full range —L < < L. 

theory. Here, x is the light-front momentum fraction of the quark and b± is the magnitude 
of the transverse relative separation coordinate. In this example, the meson eigenvalue 

n 

equation is [15| 

m = M^m. (7) 
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FIG. 5: (color online) Transverse sections of a 3-D basis function involving a 2-D harmonic oscillator 
and a longitudinal mode of Eq. [5] with box boundary conditions (wavefunction vanishes at ibL). 
The quantum numbers for this basis function are given in the caption. The basis function is shown 
for positive values of x~ and is antisymmetric with respect to = 

where the complexity of the QCD interactions among constituents is summed up in the 
addition of the effective potential f/(C), which is then modeled to enforce confinement. The 
potential in the soft wall model is U{C,) = k,'^(^ + 2k,'^{J — 1) where J is the total angular 
momentum of the hadron. Using the substitution = (^^'^R{Q, k(^ = ^/Mo^p and 

L = \m\, we arrive at the transverse 2-D harmonic oscillator wave equation whose solution 
is given in Eq. [H 



There is one additional distinction between our choice of transverse basis functions and the 
solutions of the AdS/QCD model: we adopt single-parton coordinates as the basis function 
arguments while AdS/QCD adopts a relative coordinate between the constituents. Our 
selection is natural for the applications within an external cavity that we present here and is 
most convenient for enforcing the boson and fermion statistics when dealing with arbitrary 
many partons. In future work without the external cavity, we may invoke a Lagrange 
multiplier method, analogous to the method in the NCSM and NCFC approaches , to 
separate the relative motion from the total system's motion in the transverse direction. 

The solutions of the light-front equation ([71) determine the masses of the hadrons, given 
the total internal spin S, the orbital angular momenta L of the constituents, and the index 
n, the number of nodes of the wavefunction in (. For example, if the total quark spin S is 
zero, the meson bound state spectrum follows the quadratic form = Ak^Iu + L). Thus 
the internal orbital angular momentum L and its effect on quark kinetic energy play an 
explicit role. The corresponding wavefunctions of the mesonic eigensolutions describe the 
probability distribution of the qq constituents for the different orbital and radial states. The 
separation of the constituent quark and antiquark in AdS space get larger as the orbital 
angular momentum increases. The pion with n = and L = is massless for zero quark 
mass, in agreement with general arguments based on chiral symmetry. If the total spin of the 
constituents is S* = 1, the corresponding mass formula for the orbital and radial spectrum 
of the p and u vector mesons is = AK,'^{n + L + 1/2). The states are aligned along linear 
Regge trajectories and agree well with experiment. The resulting light-front wavefunctions 
also give a good account of the hadron form factors. 

The AdS/QCD model, together with light-front holography, provides a semi classical first 
approximation to strongly coupled QCD. The BLFQ approach in this paper provides a nat- 
ural extension of the AdS/QCD light-front wavefunctions to multiquark and multi-gluonic 
Fock states, thus allowing for particle creation and absorption. By setting up and diagonal- 
izing the light-front QCD Hamiltonian on this basis, we incorporate higher order corrections 
corresponding to the full QCD theory, and we hope to gain insights into the success of the 
AdS/QCD model. 



IV. CAVITY MODE LIGHT-FRONT FIELD THEORY WITHOUT INTERAC- 
TIONS 



For a first application of tlie BLFQ approacli, we consider a non-interacting QED system 
confined to a transverse harmonic trap or cavity. For simplicity, we take the spin 1/2 leptons 
as massless. The basis functions are matched to the trap so we implement a transverse 
2-D harmonic oscillator basis with length scale fixed by the trap and finite modes in the 
longitudinal direction with APBC Since we are ultimately interested in the self-bound states 
of the system, we anticipate adoption of the NCSM method for factorizing the eigensolutions 
into simple products of intrinsic and total momentum solutions in the transverse direction . 
That is, with a suitable transverse momentum constraint such as a large positive Lagrange 
multiplier times the 2-D harmonic oscillator Hamiltonian acting on the total transverse 
coordinates, the low-lying physical solutions will all have the same expectation value of 
total transverse momentum squared. Therefore, following Ref. [l| we introduce the total 
invariant mass-squared for these low-lying physical states in terms of a Hamiltonian H 
times a dimensionless integer for the total light front momentum K 



+ P^P^ + const = P+P- = KH (8) 

where we absorb the constant into M^. The Hamiltonian H for this system is defined by the 
sum of the occupied modes i in each many-parton state with the scale set by the combined 
constant A^ = 2MqVL: 

g = 2M„p-^giH2gv2"- + l""l + l. (9) 
K ^-^ Xi 

i 

We adopt symmetry constraints and two cutoffs for our many-parton states. For symme- 
tries, we fix the total charge Z, the total azimuthal quantum number Mt, and the total spin 
projection S along the x~ direction. For cutoffs, we select the total light-front momentum, 
K, and the maximum total quanta allowed in the transverse mode of each many-parton 
state, Nmax- For the longitudinal modes, we select those with PBC from Eq. HI The chosen 
symmetries and cutoffs are expressed in terms of sums over the quantum numbers of the 
single-parton degrees of freedom contained in each many-parton state of the system in the 
following way: 



z 



(10) 
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Mt 
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^2ni + |mi| + 1 < A^, 



(14) 



where, for example, ki is the integer that defines the PBC longitudinal modes of Eq. H] 
for the i'^'^ parton. The range of the number of fermion-antifermion pairs and bosons is 
limited by the cutoffs in the modes {K and N^ax)- Since each parton carries at least one 
unit of longitudinal momentum, the basis is limited to K partons. Furthermore, since each 
parton carries at least one oscillator quanta for transverse motion, the basis is also limited 
to Nmax partons. Thus the combined limit on the number of partons is mm{K^Nmax)- In 
principle, one may elect to further truncate the many-parton basis by limiting the number 
of fermion-antifermion pairs and/or the number of bosons but we have not elected to do so 
here. 

We may refer to the quantity K as the inverse longitudinal harmonic resolution. We 
reason that as we increase K, higher longitudinal momenta states become available to the 
partons, thus allowing finer detail in the features of the longitudinal coordinate structure to 
emerge. 

In a fully interacting application, the actual choice of symmetry constraints will depend 
on those dictated by the Hamiltonian. For example, with QCD we would add color and fiavor 
attributes to the single particle states and apply additional symmetries such as requiring all 
many-parton states to be global color singlets as discussed below. Another example, which 
we adopt for the interacting QED example below, is the choice to conserve total Mj = 
Mt + S rather than conserving Mt and S separately. It is straightforward, but sometimes 
computationally challenging, to modify the symmetries in a basis function approach such as 
we adopt here. However, in order to approach the continuum limit (all cutoffs are removed) 
as closely as possible with limited computational resources, one works to implement as many 
of the known symmetries as possible. 
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FIG. 6: (color online) State density as a function of dimensionless state energy E from BLFQ for 
non- interacting QED in a trap with no net charge and for a selection of Nmax values at fixed K = 
The dimensions of the resulting matrices are presented in the legend. The states are binned in 
groups of 5 units of energy (quanta) where each parton carries energy equal to its 2-D oscillator 
quanta (2nj + |mj| + 1) divided by its light-front momentum fraction (xj = ki/K). The dashed line 
traces an exponential in the square root of energy that reasonably approximates the histogram at 
larger N^nax values. 

A. Basis space dimensions 

For our defined non-interacting cavity mode problem, we now illustrate the exponential 
rise in basis-space dimensions with increasing N^ax at fixed K, with increasing K at fixed 
Nmax and with simultaneous increase in both cutoffs. The first two situations involve a 
parton number cutoff defined by K and N^ax respectively. Only the case with simultaneous 
increase in cutoffs keeps the problem physically interesting at higher excitations since this 
is the only case with unlimited number of partons as both cutoffs go to infinity. 



In Fig. [6] we present the state density in BLFQ for massless QED in the zero couphng 
hmit for the case with no net charge Z = 0, i.e. for zero lepton number. Thus the cavity is 
populated by many-parton states consisting of fermion-antifermion pairs and photons. The 
chosen symmetries are M = and S = 0. We show results for K = Q at various values 
of Nmax spanning a range {N^ax = 8 — 23) . The states are grouped to form a histogram 
according to their energy calculated from the chosen Hamiltonian in Eq. [9] where we omit 
the constant preceding the summation for simplicity. Similarly, in Fig. [7] we present the 
state densities for Z = 3, Mt = and 5 = 1/2 at the same selected values of Nmax- 

Both Figs. [6] and [7| demonstrate the saturation of low-lying modes with increasing Nmax- 
That is, in each case, one may observe an excitation energy at which the state density reaches 
a value that no longer changes with increasing Nmax- The energy at which this saturation 
occurs, increases with Nmax- We show only the lower sections of some of the state density 
distributions but it is clear that all distributions must fall off at sufficiently high energy for 
fixed Nmax and K. 

In Fig. [8] we present the state density in BLFQ for QED in the zero coupling limit again 
for the case with no net charge Z = but with increasing K at fixed Nmax = 8. In this 
case the many-parton states at low energy continue to increase in number with increasing 
K. This is understandable from the definition of the Hamiltonian in Eq. [91 In particular, a 
typical fermion-antifermion state with each parton's light-front momentum fraction close to 
Xj = I achieves a low energy. Correspondingly, as one increases K, the population of states 
at low E grows since there are more pairs of values of Xi near | to employ for minimizing 
the energy. This reasoning easily extends to states with increasing numbers of partons so 
the net result is an increasing level density with increasing K at fixed low E and fixed Nmax- 

For the final example of state densities, we consider the case where both K and Nmax 
increase simultaneously. For simplicity, we remain with the Z = sector and take K = Nmax- 
The state densities for this example are presented in Fig. O Here, we observe trends similar 
to those shown in Fig. [8] where there is no saturation in state density at low energy. 

We take three cases depicted in Fig. [9] to illustrate the distribution of many-parton 
states over the sectors of the Fock-space. The distributions for the Nmax = -ft' = 8, 10 and 
12 examples are shown in Table [B With increasing cutoff, there is a rapid growth in the 
number of basis states within each Fock space sector. Overall, there is approximately a 
factor of 20 increase in the total many-parton basis states with each increase of 2 units in 




FIG. 7: (color online) State density as a function of dimensionless state energy E from BLFQ for 
non-interacting QED in a trap with net charge of 3 and for a selection of Nmax values at fixed 
K = Q . The dimensions of the resulting matrices are presented in the legend. The states are 
binned in groups of 5 units of energy (quanta) where each parton carries energy equal to its 2-D 
oscillator quanta (2nj + |mj| + 1) divided by its light-front momentum fraction {xi = ki/K). The 
dashed line traces an exponential in the square root of energy that reasonably approximates the 
histogram at larger N^ax values. 

the coordinated cutoff. 

Specific cases in Table [T] where no basis states may appear in a given Fock space sector 
may seem puzzling at first. However, they are understandable once the symmetries and 
constraints are examined. For example, with Nmax = K = 8 there are no states with 4 // 
pairs since the Pauli principle excludes more than 2 pairs from occupying the lowest Nmax 
and K modes. Since two // pairs must be in higher modes, either the total = 8 or 
Nmax = 8 constraint will be violated by having a total of 4 // pairs. 
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FIG. 8: (color online) State density as a function of dimensionless state energy E from BLFQ for 
non- interacting QED in a trap with no net charge and for a selection of K values at fixed Nmax = 8. 
The dimensions of the resulting matrices are presented in the legend. The states are binned in 
groups of 5 units of energy (quanta) where each parton carries energy equal to its 2-D oscillator 
quanta (2nj + |mj| + 1) divided by its light-front momentum fraction {xi = ki/K). 



All our level density results are shown as a function of the dimensionless energy. For the 
non-interacting theory in BLFQ only the kinetic term of the Hamiltonian contributes and 
the scale is available through an overall factor = 2MqQ as described above. Without 
interactions and the associated renormalization program, one cannot relate the scales at 
one set of {K, Nmax) values to another. Ultimately, one expects saturation will arise with 
interaction/renormalization physics included as one increases the set of {K, N^ax) values. 

These state densities could serve as input to model the statistical mechanics of the system 
treated in the microcanonical ensemble. Of course, interactions must be added to make the 
model realistic at low temperatures where correlations are important. After turning on 
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FIG. 9: (color online) State density as a function of dimensionless state energy E from BLFQ for 



non-interacting QED in a trap with no net charge and for K 



The dimensions of the 



resulting matrices are presented in the legend. The states are binned in groups of 5 units of energy 
(quanta) where each parton carries energy equal to its 2-D oscillator quanta (2nj + |mj| + 1) divided 
by its light-front momentum fraction (xj = ki/K). 

the interactions, the challenge will be to evaluate observables and demonstrate convergence 
with respect to the cutoffs {N^ax and K). Independence of the basis scale, f2, must also 
be obtained. These are the standard challenges of taking the continuum limit. We will 
address these topics in a separate investigation. For the current effort, we present a smooth 
representation for selected histograms, an exponential fit adopted from the well-known Bethe 
formula. 



p(^) = 6exp(Va^), (15) 
where the precise values of the fitted constants are provided in the legends. We provide 



// pairs / bosons 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Total 











210 





1122 





67 





1 








1400 










495 





11318 





2936 





69 





1 


14819 










1001 





73600 





63315 





4027 





69 


142013 


1 


420 


1932 


8190 


1040 


588 


8 


2 














12180 




990 


10512 


86856 


33632 


36672 


1604 


640 


8 


2 








170916 




2002 


40810 


574860 


503040 


929064 


99962 


60518 


1770 


644 


8 


2 


2212680 


2 


5961 


1560 


1133 


4 


1 




















8659 




64240 


59240 


97584 


4040 


1513 


4 


1 














226622 




427730 


942240 


2806624 


381608 


249825 


4928 


1565 


4 


1 








4814525 


3 


218 
































218 




25584 


1528 


554 


























27666 




808034 


222336 


200676 


2592 


602 




















1234240 


4 








































16 
































16 




19325 


168 


20 


























19513 



TABLE I: Number of many-parton basis states in each Fock-space sector for two of the Nmax = K 
cases depicted in Fig. [9l The counts are organized according to the number of fermion-antifermion 
(//) pairs and the number of bosons in each sector. The first hne in each // row corresponds 
to the Nmax = K = 8 case which has a total of 22,457 states, while the second line corresponds 
to the Nmax = K = 10 case which has a total of 440,039 states. The third line in each // row 
corresponds to the N^ax = K = 12 case which has a total of 8,422,971 states. In this last case, 
there is a single 12-boson state not listed to save space. The last column provides the total for 
that row. 



these exponential fits in Figs. [6] and [7] where the low- lying state density exhibits saturation. 
B. Distribution functions 

In order to illustrate the potential value of the BLFQ approach, we present light-front 
momentum distribution functions for two simple toy models, based on results presented 
in Fig. [9l In the first example, we consider a model for a weak coupling regime and, in 
the second example, we consider a model for strong coupling behavior. In both cases we 
introduce a simple state that is an equally-weighted superposition of basis states. In the 
weak coupling case, we retain all basis states below a cutoff {Ecut = 25) in the dimensionless 
energy scale of Fig. [9] for a given value of = N^ax- That is, we imagine a situation 
where only the low-lying unperturbed many-parton basis states mix equally to describe a 
low-lying physical state of a weakly-coupled physical system. In the strong coupling case we 
retain all basis states of Fig. [9] with equal weights for a given value oi K = Nmax- Here, we 
imagine the coupling is so strong as to overwhelm the unperturbed spectrum and to produce 
a simple low-lying physical state with equal admixtures of all available basis states. These 
states, labeled I^E^^,) and {"^s), where the w ("s") represents "weak" ("strong") respectively, 
are written as normalized sums over their respective sets of many-parton basis states 
as 



where "a" represents "w" or "s" and the sum runs over the Da respective many-parton 
states. For our present application to probability distribution functions, the phases of the 
individual terms in expansion are not relevant so we choose all of them to be positive for 
simplicity. 

Selected light-front momentum distributions n{x) for these two model states are shown 
in Figs. [To] and [m The fermion and antifermion distributions are the same in these limiting 
examples. Light-front momentum distributions are probability distributions emerging after 
integration over transverse degrees of freedom. With our present selection of basis states, 
the light-front momenta take discrete values leading to discrete-valued distributions (his- 
tograms). However, for convenience, we present smooth distributions in Figs. [TOl and [TT] 




(16) 



generated by spline interpolations. 

The parton distributions at fixed Nmax = K satisfy both the normalization condition: 



The index i runs over the parton species (fermion, antifermion, boson) and the index k runs 
over the discrete values of light-front momenta corresponding to the integers in Eq. H] where 



The top panels of Figs. [lO] and [TT] display the light-front momentum distributions at 
Nmax = -ft' = 8 for the "weak" and "strong" coupling models, respectively. The lower panels 
present the boson distribution functions for three Nmax = K values ranging from 8 to 12 for 
the same models. 

The fermion distributions are found to track the boson distributions with increasing 
Nmax = K and are not shown in the lower panels. We also comment that the total mo- 
mentum distribution fractions carried by the separate parton species appear approximately 
independent of Nmax = K over the range 8 — 12. About two-thirds of the total light-front 
momentum is carried by the fermions plus antifermions. This division is characteristic of 
both the weak and strong coupling models over the Nmax = -ft' = 8 — 12 range we examined. 

The top panel of Fig. [11] indicates a peak in the vicinity of the minimum light-front 
momentum fraction carried by a single parton in this basis, x = |, for both the fermions 
and the bosons. This appears to be a characteristic of this strong coupling toy model and 
is illustrated in the lower panel of the same figure where the peaks in the boson light-front 
momentum distributions appear to track well with the inverse of Nmax = K. Clearly, with 
this toy model the distribution functions do not converge with increasing Nmax and K. 

For comparison, we note that with the weak coupling toy model the peaks of the boson 
distributions shown in Fig. [TD] appear to be stable with increasing Nmax and K, and the 
distribution function appears to be reasonably well converged at Nmax = K = 12. Based on 
these observations we anticipate good convergence for weakly interacting theories like QED. 

The lack of convergence of our strong coupling toy model may be worrisome for applica- 
tions in QCD, but one should keep in mind that this toy model is far from realistic: all basis 
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and total light-front momentum conservation 
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FIG. 10: (color online) Light front momentum distribution functions for states representing a 
weak coupling paradigm. The top panel displays the distributions at Nmax = K = 8. The 
antifermion distribution is the same as the fermion distribution. The total momentum fraction 
carried by the fermion plus antifermion distribution is 0.66 while the boson distribution carries 
the remaining fraction 0.34. The bottom panel displays the boson distributions at three different 
values of Nmax = K that are labeled. 





FIG. 11: (color online) Light front momentum distribution functions for states representing a 
strong coupling paradigm. The top panel displays the distributions at Nmax = K = 8. The 
antifermion distribution is the same as the fermion distribution. The total momentum fraction 
carried by the fermion plus antifermion distribution is 0.65 while the boson distribution carries 
the remaining fraction 0.35. The bottom panel displays the boson distributions at three different 
values of Nmax = K that are labeled. 



states are retained with equal weight. Nevertheless, it is interesting to consider the trends of 
this model with increasing Nm.ax = K. For background, one may recall that deep-inelastic 
lepton scattering from a hadron in the scaling region Q"^ oo provides a measure of the 
hadron's charged quark distribution functions. With more detailed resolution provided by 
the virtual photon exchange (increasing Q leads to shorter wavelengths), experiments re- 
veal that the charged quark distributions evolve to lower values of light-front momentum 
fraction, x. The pattern shown in the lower panel of Fig. [11] with increasing N^ax = K 
is reminiscent of this experimental trend with increasing Q. Given the simplicity of the 
strong interaction model, one may infer that the evolution of multi-parton phase space with 
increasing N^ax = K could play a significant role in the evolution of light-front momentum 
distribution functions with improved resolution through increasing Q. 



C. Extension to color 



We can extend the approach to QCD by implementing the SU(3) color degree of freedom 
for each parton - 3 colors for each fermion and 8 for each boson. For simplicity, we restrict 
the present discussion to the situation where identical fermions occupy distinct space-spin 
single-particle modes. The case where we allow multiple space-spin occupancies by identical 
fermions leads to color space restrictions. We will address this additional complexity in a 
subsequent investigation. 

We consider two versions of implementing the global color-singlet constraint for the re- 
stricted situation under discussion here. In both cases we enumerate the color space states 
to integrate with each space-spin state of the corresponding partonic character. 

In the first case, we follow Ref. 18| by enumerating parton states with all possible values 
of SU(3) color. Thus each space-spin fermion state goes over to three space-spin-color states. 
Similarly, each space-spin boson state generates a multiplicity of eight states when SU(3) 
color is included. We then construct all many-parton states having zero color projection. 
Within this basis one will have both global color singlet and color non-singlet states. The 
global color-singlet states are then isolated by adding a Lagrange multiplier term in many- 
parton color space to the Hamiltonian so that the unphysical color non-singlet states are 
pushed higher in the spectrum away from the physical color single states. To evaluate 
the increase in basis space dimension arising from this treatment of color, we enumerate 
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FIG. 12: (color online) Number of color space states that apply to each space-spin configuration 
of selected multi-parton states for two methods of enumerating the color basis states. The upper 
curves are counts of all color configurations with zero color projection. The lower curves are counts 
of global color singlets. 



the resulting color-singlet projected color space states and display the results as the upper 
curves in Fig. [121 

In the second case, we restrict the basis space to global color singlets and this results in 
the lower curves in Fig. [T^l The second method produces a typical factor of 30-40 lower 
multiplicity at the upper ends of these curves at the cost of increased computation time for 
matrix elements of the interacting Hamiltonian. That is, each interacting matrix element in 
the global color-singlet basis is a transformation of a submatrix in the zero color projection 
basis. Either implementation dramatically increases the state density over the case of QED, 
but the use of a global color-singlet constraint is clearly more effective in minimizing the 
explosion in basis space states. 



We note that, for the pure multi-fermion basis space sector, shown in the upper left panel 
of Fig. [121 "we could have produced the lower curve using methods introduced and applied 



successfully in 1 + 1 dimensional QCD [19|. That is, the number of global color singlets 
for a given fermion-only basis state, with other (non-color) quantum numbers specified, is 
independent of the number of spatial dimensions provided there is at least one. 

V. CAVITY MODE LIGHT-FRONT FIELD THEORY - ELEMENTS OF THE IN- 
TERACTING THEORY 

Now we briefly address the interacting theory where a primary concern will be to manage 
the divergent structure of the theory. There are two possible locations for divergences in 
a Hamiltonian basis function approach: (1) the matrix elements themselves diverge, or (2) 
the eigenvalues diverge as one or more cutoffs are removed. 

In our cavity field theory applications with interactions, we manage these divergences 
with the help of suitable counterterms, coupling constant and mass renormalizations, and 
boundary condition selections. The development of counterterms is expected to be straight- 
forward as seen, for example, in Ref. The infrared divergences in light-front momentum 
arising in both the fermion to fermion-boson vertex as well as in the instantaneous fermion- 
boson interaction are expected to be well-managed by previously defined counterterms 
suitably transcribed for the transverse basis functions we have adopted. We anticipate this 
prescription will work since the longitudinal modes we adopt are similar to those used in 
Ref. 3|. Finally, we expect the transverse ultraviolet divergences to be suitably-managed 
with our basis function selection. 

Since we are introducing a basis-function approach for the transverse degrees of freedom, 
we need to investigate convergence rates with increased cutoff of the transverse modes N^ax- 
Here, in a simple set of examples, we outline how we can search for additional sources of 
divergence with the help of a perturbation theory analysis. 

Consider the second order energy shift, AE, induced on a single parton in the transverse 
mode {n,m) by its coupling V to partons in higher energy transverse modes {n',m'): 




(19) 
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where we use the following notation and properties of the 2-D harmonic oscillator 
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n' = 2n' + \m'\ 
= {n' + 1)0 
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and we converted the sum to an integral taking the degeneracy into account with p. 

Thus, according to perturbation theory, we expect a UV divergence if the matrix element 
falls off too slowly with increasing n'. In particular, if the matrix element falls approximately 
as (n')~^, then we expect a logarithmic divergence since the integrand will have a net {n')~^ 
dependence. If the falloff is even slower then we encounter a more serious divergence. 

Another possible source of a log divergence could arise within the selected sum over m' in 
which case p, the level density factor in the integrand, is unity. Then, if the matrix elements 
for fixed n, n' are approximately constant with increasing m', we again find a log divergence 
in the sum over m' . 

For a first investigation, we have examined the behavior of various sets of matrix ele- 
ments for the fermion to fermion-boson vertex. For the purpose of this investigation, we 
adopt periodic (antiperiodic) boundary conditions for the longitudinal modes of the bosons 
(fermions) and we hold spin projections fixed for initial and final states. We then adopt spe- 
cific values for the longitudinal momentum fractions observing the conservation rule. The 
trends we examine should not be sensitive to the specific values adopted. 

To date, we have found no matrix element trends with increasing transverse energy that 
would imply new divergences. All such matrix elements sets we examined, within a per- 
turbative analysis, fall faster than the inverse square root of the principal quantum number 
n' . Also, we found no sets that remained constant with increasing m' (and thus m) holding 
n, n' fixed. 

We portray in Figs. [13] and [14] representative sequences of how these off-diagonal matrix 
elements behave as one increases the difference in the initial and final state principal quan- 
tum numbers. We also portray two interesting cases where the fermion and fermion-boson 
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FIG. 13: (color online) Behavior of representative fermion to fermion-boson matrix elements in 
BLFQ. The quantum numbers specifying the parton transverse modes (nj,mj) in the matrix el- 
ements are given in the legend. Only the transverse mode contributions to the matrix elements 



are shown. Results are also shown with a multiplicative factor of y/n + 1 applied to help search 
for a logarithmic divergence by obtaining a resulting flat behavior, when it occurs. Overall matrix 
element normalization depends on the specific values of light-front momentum fractions carried by 
the interacting partons. 

principal quantum numbers track each other, cases that do not enter a perturbative analy- 
sis. Note that we have limited the illustrations to the transverse components of our matrix 
elements. We also select cases where the fermion spin is flipped, cases that are proportional 
to the fermion mass. We set the fermion mass to unity so the results are expressed in units 
of the fermion mass. 

We further limit our presentation to the case where all partons remain in the orbital 
projection quantum number zero state and the 2-parton (fermion-boson) states have each 
parton in the same transverse state. For the non-perturbative illustrative cases, we display 
the matrix element trends where all partons remain in the same transverse mode, (n, 0). 
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FIG. 14: (color online) Behavior of representative fermion to fermion-boson matrix elements in 
BLFQ. The quantum numbers specifying the parton transverse modes (nj,mj) in the matrix ele- 
ments are given in the legend. Only the transverse mode contributions to the matrix elements are 



shown. Results are also shown with a multiplicative factor of y/n + I applied to identify logarithmic 
divergence by a resulting flat behavior, when it occurs. Note that for one of the cases shown, the 
resulting matrix elements vanish with increasing n while the other case shows a constant trend but 
does not enter a second-order perturbative analysis. Overall matrix element normalization depends 
on the specific values of light-front momentum fractions carried by the interacting partons. 

When the single fermion state has an even value of the principal quantum number n as 
shown in Fig. [131 the matrix elements appear to be well-behaved either when the 2-parton 
configuration is the lowest accessible case ((0000|) or when the each of the two partons 
resides in the state with the same principal quantum number as the single fermion state. 
We demonstrate anticipated good convergence with increasing n by showing that the matrix 



elements, when multiplied by \/n + 1, still fall with increasing n. 

For the case when the single fermion state has an odd value of the principal quantum 
number n as shown in Fig. [T31 the situation is somewhat different. For the matrix element 



set entering a perturbative analysis, the matrix elements fall to zero with increasing n 



sufficiently fast that multiplying by + 1 does not significantly distort the trend to zero. 
However, the large n behavior of the fermion-boson matrix element, with all partons at the 
same n, is seen to go approximately as \/n + 1. This is best seen in Fig. [H] where the 
matrix elements are multiplied by + 1 and the result appears to be a nonzero constant 
at large n. Since this trend does not appear in a second order perturbation theory analysis, 
we must await the full Hamitonian diagonalization in sufficiently large basis spaces to better 
understand its role in the convergence with increasing Nmax- 

As a result of this initial analysis, we anticipate that straightforward adoption of coun- 
terterm methods previously introduced will be sufficient for managing the identified 
divergences in BLFQ. However, it seems advisable to have an alternative scheme for com- 
parison. Therefore, we plan to adopt a second approach that involves a recently proposed 
sector-dependent coupling constant renormalization scheme {2^. Another alternative which 
we may adopt uses the Pauli-Villars regulator 211 ]. 

Even without regulating the possible divergencies in the Hamiltonian, it is possible to 
get some idea how the cutoffs for the basis space dimensions, as discussed in the previous 
section, affect the eigenvalue spectra of the Hamiltonian. In Fig. [15] we show the eigenvalues 
(multiplied by K) for a light-front QED Hamiltonian in a basis space limited to the fermion 
and fermion-boson sectors. For this particular example we chose the harmonic oscillator 
parameters as = 0.1 MeV and Mq = 0.511 MeV, and the fermion mass was chosen to be 
equal to Mq. The interaction terms include the fermion to fermion-boson vertex and the 
instantaneous fermion-boson interaction. We chose the basis space such that the basis states 
have total Mj = Mt + S* = |, and we simultaneously increase the K and Nmax cutoffs. As 
a result, the size of the Hamiltonian matrices increase rapidly. For K = Nmax = 2, 3, 4, 5, 6, 
the dimensions of the corresponding matrices are 2x2, 12 x 12, 38 x 38, 99 x 99 and 208 x 208 
respectively. 

The number of the single fermion basis states increase slowly with increasing K = Nmax 
cutoff. For K = Nmax = 2, 3, 4, 5, 6 the number of single fermion basis states is 1, 2, 2, 3, 3, 
respectively. Our lowest-lying eigenvalues correspond to solutions dominated by these states 
and they appear with nearly harmonic separations in Fig. [15] as would be expected at the 
coupling of QED. 

The higher eigenstates are the ones dominated by the fermion-boson basis states that in- 



teract with each other in leading order through the instantaneous fermion-boson interaction. 
Their multiphcity increases rapidly with increasing K = Nmax and they exhibit significant 
mixing with each other as well as weak mixing with the lowest-lying states. The eigenvalues 
dominated by the fermion-boson basis states cluster in nearly degenerate groups above the 
lowest-lying states. 

Further progress requires that we implement our renormalization program as well as a 
major expansion in the basis space size. These efforts are underway and will be presented 
in a subsequent paper. 

VI. CONCLUSIONS AND OUTLOOK 

Following successful methods of ab initio nuclear many-body theory, we have introduced 
a basis light-front quantization (BLFQ) approach to Hamiltonian quantum field theory and 
illustrated some of its key features with a cavity mode treatment of massless non-interacting 
QED. 

Cavity mode QED, with a 2D harmonic oscillator for the transverse modes and longitudi- 
nal modes chosen with periodic boundary conditions, exhibits the expected dramatic rise in 
many-parton basis states as the cutoffs are elevated. With the non-intracting cavity-mode 
Hamiltonian, we obtain the state density distributions at various choices of the regula- 
tors. These basis state densities provide initial elements of a quantum statistical mechanics 
approach to systems treated in the BLFQ approach. We then illustrated the access to light- 
front momentum distribution functions in this approach with simple models of wavefunctions 
that reflect possible interaction effects. 

In order to extend our method to QCD, we have evaluated two methods for treating the 
color degree of freedom. Since large sparse matrices will emerge, we argue that it is more 
efficient in storage requirements to adopt multi-parton basis states that are global color 
singlets and we presented sample measures of the efficiency gains over basis states with 
color-singlet projection alone. To achieve this savings in storage (reduced matrix size) we 
will incur an increase in the computational effort for the non-vanishing matrix elements. 

We have also outlined our approach to managing the expected divergences that will pre- 
serve all the symmetries of the theory. An initial inspection of the interaction vertices of 
QED in the BLFQ approach shows smooth behaviors that, following a second-order pertur- 
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FIG. 15: (color online) Eigenvalues (multiplied by K) for a non-renormalized light-front QED 
Hamiltonian which includes the fermion-boson vertex and the instantaneous fermion-boson inter- 
action without counterterms. The basis is limited to fermion and fermion-boson states satisfying 
the symmetries. The cutoffs for the basis space dimensions are selected such that K increases 
simultaneously with the N^ax- The harmonic oscillator parameters were chosen as Vt = 0.1 MeV 
and Mq = 0.511 MeV, and the fermion mass was chosen to be equal to Mq. 



bative analysis, are not expected to lead to divergences. It appears that the cavity mode 
treatment, with the type of basis spaces we have selected, will encounter the divergences 
in a more subtle fashion as cutoffs are elevated. For illustration, we presented an initial 
example in QED, without renormalization, of the mass spectrum for a single electron in a 
transverse cavity coupled to single photon modes. The computational requirements of the 
BLFQ approach are substantial, and we foresee extensive use of leadership-class computers 
to obtain practical results. 
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